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Properties of atomic nuclei important for the prediction of astrophysical reaction rates 
are reviewed. In the first part, a recent simulation of evolution and nucleosynthesis of 
stars between 15 and 25 M is presented. This study is used to illustrate the required 
nuclear input as well as to give examples of the sensitivity to certain rates. The second 
part focusses on the prediction of nuclear rates in the statistical model (Hauser-Feshbach). 
Some of the important ingredients are addressed. Discussed in more detail are approaches 
to predict level densities, parity distributions, and optical a-nucleus potentials. 

1. Introduction 

The knowledge of nuclear reactions is crucial to model the evolution of stars and to 
determine the amounts of produced nuclei. From the abundance patterns of the elements 
we can learn something about how they were created and consequently it is possible to 
study the conditions in supposed astrophysical sites and, finally, to attempt to trace the 
origin and history of the Universe. An in-depth knowledge of the sub-atomic processes is 
required in order to understand many large-scale effects. 

The demands of astrophysics challenge our ability to describe and predict nuclear prop- 
erties and nuclear reactions. Due to high temperatures and densities in the stellar plasma, 
a large number of unstable nuclides are involved in the nucleosynthetic processes. Some of 
their radioactive decay lines can be observed by modern satellite observatories and provide 
an important tool to test the hydrodynamics of stellar models by comparison. However, 
this task can only be performed if the producing nuclear reactions are sufficiently well 
known. In this respect, not only reactions involving unstable targets are important but 
also such along the line of stability. Despite of temperatures in the range of millions to bil- 
lions of Kelvin, the respective energies of the interacting particles are quite low by nuclear 
physics standards, from thermal energies up to a few MeV. Due to the size of the reaction 
networks (the nucleosynthesis calculations for massive stars, which are described below, 
include about 2400 nuclides in more than 15000 reactions), theoretical predictions will al- 
ways remain important. However, the relevance of the low-energy region poses problems 
to both experimental and theoretical approaches. For neutron-induced reactions, to re- 
solve the transition between and interplay of different reaction mechanisms - from direct 
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to resonant interactions and finally to the Hauser-Feshbach regime - remains a challenge. 
Except for activation experiments, standard techniques only measure the resonant part 
of the cross section. On the other hand, energy and strength of resonances are difficult 
to predict with the required accuracy. The situation is somewhat improved whenever the 
statistical model of nuclear reactions (Hauser-Feshbach theory) can be applied, which uses 
resonance averages. Fortunately, the majority of reactions in astrophysics involving the 
strong interaction can be described in this approach [UT]. When predicting cross sections 
for astrophysical applications in such a way, slightly different points are emphasized than 
in pure nuclear physics investigations. Since most of the ingredients for the calculations 
are experimentally undetermined (in some cases even for nuclides at or close to stabil- 
ity), one has to develop reliable phenomenological or microscopic models to predict these 
properties with an acceptable accuracy across the nuclear chart. Therein one has to be 
satisfied with a more limited accuracy as compared to usual nuclear physics standards. 
Considering the substantially larger uncertainties in many astrophysical scenarios, this 
seems to be adequate. For certain selected reactions, however, the sensitivity of the as- 
trophysical results is so high that an accuracy of 10% or better is necessary. Examples of 
this are the reactions 12 C(a,7) 16 and (a,n) and (0,7) on 22 Ne, as discussed below. When 
studying details in s-process branchings, cross sections have to be determined to 1%. This 
is currently only possible with experiments or by a combination of experiment and theory, 
e.g. when supplementing resonance measurements with direct capture calculations. 

In this review I first present new calculations of the evolution and nucleosynthesis 
of massive stars, involving reaction networks of unprecedented size. A non-exhaustive 
selection of important nuclear reactions is then given and their importance is shown 
by the impact on the results. In the second part of the paper, the current standard 
approach to predict nuclear cross sections and reaction rates is outlined and the remaining 
uncertainties and challenges are illustrated by a few examples. 

2. Nucleosynthesis in massive stars 
2.1. The model 

Stars above ~ 10 M & are responsible for producing most of the oxygen and heavier 
elements found in nature. Numerous studies have been devoted to the evolution of such 
stars and their nucleosynthetic yields. However, our knowledge of both the input data 
and the physical processes affecting the evolution of these stars has improved dramatically 
in recent years. Thus, it became worthwhile to improve on and considerably extend the 
previous investigations of pre- and post-collapse evolution and nucleosynthesis. The first 
calculation to determine, self-consistently, the complete synthesis of all stable nuclides 
in any model for a massive star [SB] is discussed here. The calculations were performed 
using the stellar evolution code KEPLER [TH] with several modifications relative to [TH] 
(mass loss due to stellar winds, improved adaptive network) and updates (OPAL95 opacity 
tables, neutrino loss rates). According to the topic of this report, the focus is on giving an 
outline of the updates concerning the nuclear reactions involving the strong interaction. 
For further details of the calculations and results the reader is referred to the full paper 

As in jTH], two reaction networks are used. A small network directly coupled to the 




Figure 1. Postexplosive production factors (vertical axis) versus mass number (horizontal 
axis) of a population I star with 15 M progenitor mass. The horizontal axis Comparison 
is relative to solar abundances. The relative overproduction of 61>62 Ni may indicate some 
inherent nuclear uncertainty in the (11,7) destruction cross section. See text and the full 
paper [HE] for details. 



stellar model calculation provides nuclear energy generation, i.e. it is solved implicitly for 
each time-step in each zone. This smaller network is essentially the same as in |7B] . but 
with updated nuclear rates as described in the following. 

One of the major improvements over [TH] and other stellar models is that, for the first 
time, the synthesis of all nuclides of any appreciable abundance is followed simultaneously 
in an adaptive network of unprecedented size. Using a library containing rate information 
for 4,679 isotopes from hydrogen to astatine, the "adaptive" network automatically adjusts 
its size to accommodate the current nuclear flows. This saves CPU time and thus allows 
to perform the calculations within reasonable time. Because of convective coupling of 
zones, the same network must be used throughout the star. 

Within such a network, the impact of nuclear rates on burning and evolution of the 
star can be studied fully self-consistently. 

2.2. Nuclear input 

The most extensive published library of theoretical reaction rates to date |59|60j was 
used as the backbone of the reaction rate sets. Details of the predictions are given in Sec. 
01 For the network described here, the rates based on the FRDM set were implemented. 
This theoretical set was supplemented with experimental neutron capture rates along the 
line of stability b\. Experimental (0,7) rates were implemented for 70 Ge [TH] and 144 Sm 
[BTj . The derived a+ 70 Ge and «+ 144 Sm potentials were also utilized to recalculate the 
transfer reactions involving these potentials. Semi-empirical rates were implemented for 
a-capture reactions on self-conjugated (N = Z) nuclides (HSI- For the important rate 
12 C(a,7) 16 we used an updated rate (5(300) = 146 keV barn = 1.2x Ref. QUI) and 
temperature dependence ^UJ- Similarly important is the (a,n)/(a,7) branching [HI] on 
22 Ne (see |ZS])- Table 1 defines the standard rate set. For comparison, additional sets of 
experimental and theoretical rates were used for elements below neon: Refs. |76|28j and 
NACRE [3 . Experimental @~, j3 + , and a-decay rates were taken from [70] and theoretical 
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Table 1 

Nuclear reaction rate inputs (for details see Tables 1 and 2 in Ref. 



Reference 



Type 



Bao et al. (2000) [5J 

Buchmann (1996) [Til] 

Caughlan & Fowler (1988) [JJ] 

Fowler, Caughlan & Zimmerman (1975) [iT| 

Fiilop et al. (1996) [U] 

Giessen et al. (1994) [23] 

Harris et al. (1983) [2l] 

Hansper et al. (1989) |23] 

Iliadis et al. (2001) [3T] 

Kappeler et al. (1994) [54] 

Kiener et al. (1993) [35] 

Landre et al. (1990) |3J] 

Mitchell et al. (1985) |35| 

Morton et al. (1992) [517] 

Rauscher & Thielemann (2000,2001) |5SI50| 

Rauscher et al. (1994) 

Rauscher et al. (2000) 

Somorjai et al. (1998) [HZ] 

Scott et al. (1992) [55] 

Sevior et al. (1986) [02] 

Scott et al. (1991) [H3] 

Wrean et al. (1994) [77] 

Wiescher & Kettner (1982) [73] 

Winters & Macklin (1988) [74] 



0,7) 

12 C(a,7), modified (see Woosley, this volume) 

light targets, charged projectiles 

light targets, charged projectiles 

™Ge(a, 7 ) 

18 0(a, 7 ) 

light targets, charged projectiles 

45 Sc(a,p), 45 Sc(a,n) 

(P,7) 

22 Ne+a, modified (see Woosley, this volume) 
13 N(p, 7 ) 

14 N(a,p), 17 0(p, 7 ) 
42 Ca+a 
48 Ti(a,p) 
theory 

light targets, neutron- and a-induced 

a-capture on isospin symmetric targets 

144 Sm(a,7) 

34 S+a 

38 Ar+a 

41 K(a,p) 

9 Be(a,n) 

15 0(a,p) 

20 Ne(n, 7 ) 



(3 and j3 + rates from [ISj- As a special temperature-dependent 180 Ta decay [E] 

was implemented. The production ratio 180m Ta/ 180 Ta had to be computed offline (after 
production and ejection) because the isomer was not included in the network as a separate 
species. The estimate of this ratio is based on the data of jB] and a derivation given in 
Appendix B of [55] , 

For A < 40 recent theoretical weak rates [H] were also included. The i/-process was 
not followed for nuclides with Z or iV larger than 40 and neither was a possible r-process 
in high-entropy layers close to the surface of the proto-neutron star. However, a slight 
n-process could be found, due to high neutron flux generated at the base of the He-shell 

123 ■ 

2.3. A non-exhaustive selection of important reactions 

It should be noted that only a selection of a few important rates is presented in this 
section. If a reaction does not show up here, this does not imply that it is of no importance. 
It just means that there was not sufficient space to discuss the reaction properly. Main 
reactions in the quiescent burning phases of a star (e.g. 12 C+ 12 C in carbon burning) as 
well as those many neutron capture reactions necessary for a detailed study of s-process 
branchings will always remain on the list of important reactions. However, a few reactions 
were chosen subjectively for discussion here. 
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2.3.1. 12 C(a, 7 ) 16 

It is well known that the leading nuclear uncertainty afflicting modern studies of stellar 
evolution and nucleosynthesis continues to be the reaction rate of 12 C(a,7) 16 0. This re- 
action competes for the consumption of a-particles with the triple-alpha reaction during 
helium burning. It determines the ratio of C to O at the onset of the subsequent burning 
stage, i.e. carbon burning. This has important implications not only for nucleosynthesis 
but determines the further evolution of the star and even the nature of the collapse and 
explosion. 

Weaver & Woosley [72] suggested a preferred value of 170±20 keV b for the total 
S-factor at 300 keV, based on nucleosynthesis arguments. Recently, those calculations 
were repeated using more modern stellar models ,26 and arrived at basically the same 
conclusions as [72]. For that reason, the currently preferred choice for the absolute value 
and energy dependence of the rate is the one given in Sec. 12.21 

A further implication of the nucleosynthesis study is that the rate has to be known to an 
accuracy of < 10%. Recent measurements |10|11|9|38I6 9] give recommended values in the 
range 145 to 165 keV b but with considerable uncertainty. Only the work of Tischhauser 
et al. [HH] is beginning to approach that precision. 

2.3.2. 22 Ne(a,7) 26 Mg and 22 Ne(a,n) 25 Mg 

Behind 12 C(a,7) 16 0, the 22 Ne(a,n) rate and the (a,n)/(a,7) branching at 22 Ne become 
the second largest nuclear uncertainties in calculating the nucleosynthesis in massive stars. 
The weak s-process component, including a large number of nuclides produced by neutron 
captures up to about A = 90, is made in massive stars. The relative production of these 
nuclei compared to other abundant species like O and to one another is sensitive to the 
cross sections of neutron poisons like 25 Mg and 16 [73], and to the rate of 22 Ne(a,n) 25 Mg. 

The (a,n) rate has always been quite uncertain [3|. Recent experimental work |3*2~| 
exhibits higher accuracy but is still uncertain enough to accommodate a factor of two 
uncertainty in many important s-process products for 80 < A < 90 (see Fig. |2J). 

2.3.3. The 62 Ni(n,7) case 

This case is a good example of the difficulties encountered when trying to predict 
reaction rates for final nuclei with low level densities. For neutron-induced reactions at 
low energies, close to magic numbers, and far off stability where low separation energies are 
encountered, a problem emerges. In such targets, the level density is too low to allow the 
application of the statistical model [HI]. Also for other nuclides it is not straightforward 
to bridge the region of thermal energies to the region of overlapping resonances where the 
Hauser-Feshbach formalism can be used. Single resonances and direct reactions become 
important. This is also an issue for neutron-rich nuclei in the r-process path with low 
neutron-separation energies. 

With the reaction 62 Ni(n,7) 63 Ni it was previously attempted to extrapolate thermal 
data to s-process energies of up to a few hundred keV. Two compilations give disagreeing 
30 keV cross sections |4l5j . based on the same thermal data. Both extrapolations assume 
s-wave behavior of a direct capture component. The more recent one includes a sub- 
threshold resonance contributing to the thermal cross section. 

A calculation of the direct capture component using DWBA found a considerable p- 
wave contribution which enhances the cross section at 30 keV [57 . Thus, even when 
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Figure 2. Nucleosynthesis of the weak s-process with varied rate of 22 Ne(a,n) 25 Mg in a 
25 M Q star of initial solar metallicity. A factor of two overproduction (dotted lines) is 
acceptable because lower mass stars make less s-process and low metallicity stars make 
very little. 



including the subthreshold resonance, the 30 keV value is coincidentally similar to the 
value in the older compilation (Fig. |2J). However, also the general energy dependence of 
the cross section is altered. Resonances were also included but they only contribute less 
than 15%. The enhanced cross section has an important impact on s-processing in massive 
stars. The previously seen overproduction of 62 Ni in stellar models (see Fig. ^) can be 
cured when using the enhanced rate because of increased destruction of this nucleus with 
the larger neutron capture rate |5b|57j . An experimental verification of this result would 
be desireable. 

Further neutron capture reactions on Ni isotopes might also be suspect but they do not 
show the additional complication of a sub-threshold resonance. 
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Figure 3. Direct neutron capture Maxwellian averaged cross section of 62 Ni. The final 
value is given by adding the resonant contribution to the "total" direct term. Upper and 
lower limits on the direct components are from experimental errors on the input, i.e. in 
the thermal scattering length and the spectroscopic factors. 



3. Statistical model calculations 

3.1. A reaction rate library 

A recently published large-scale reaction rate library includes neutron-, proton-, and 
a-induced reactions on all target nuclei from Ne up to Bi from proton-dripline to neutron- 
dripline [59 60J. Due to the fact that many very short-lived nuclides can be produced in 
astrophysical sites, it is necessary to provide cross sections and rates for about 4600 
targets and 32000 reactions. These numbers show that theory will always play a major 
role in providing cross sections, despite the potential of future Rare Isotope Accelerators. 
This rate library has already been adopted as a standard for calculating nucleosynthesis 
in stellar evolution and in type II supernovae |56|25j . Fits to the astrophysical reaction 
rates - ready for direct astrophysical application - as well as tables of cross sections, 
reaction rates, and nuclear inputs for all possible reactions with light projectiles can also 
be downloaded from \http://nucastro. org/reaclib.htm\ . 

3.2. Global calculations 

The calculations for the above library were performed with the Hauser-Feshbach code 
NON-SMOKER [SH] which is especially tuned to such large-scale predictions. The theoret- 
ical determination of nuclear properties in this code is based on microscopic or macroscopic- 
microscopic approaches suitable for an application far off stability. As mentioned before, 
due to the large number of unstable nuclei to be implemented in nuclear reaction networks 
for astrophysics, nuclear properties needed for the calculation of reaction rates have to 
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be predicted. Necessarily, there is a trade-off with accuracy because one cannot expect 
to predict those properties with similar accuracy as measured data. However, whenever 
measured properties can be used, this will, of course, lead to an improvement in the cal- 
culated rates. Unfortunately, this is not possible for the vast majority of reactions off the 
line of stability (and sometimes not even at stability). 

Occasionally, a misconception about the best models for nuclear properties seems to 
arise. The question is not merely one of microscopic versus phenomenological approaches, 
it should rather be "what can we learn from the different models and how do we imple- 
ment the findings in a simple way, suited for large-scale computations T . On one hand, 
first-principles microscopic models can be more satisfying from a philosophical point of 
view. Such a model will of course be most preferrable as soon as the perfect one has 
been identified. Unfortunately, the current state-of-art sees competition of many different 
microscopic approaches, each with its own advantages and drawbacks, but we are still 
far from a complete, unified picture. Thus, one has to be careful when using the term 
"microscopic" in the sense of a general remedy. 

On the other hand, the problem inherent in all microscopic approches is that they 
are CPU expensive. Although further progress regarding CPU speed and perhaps al- 
gorithmical effectiveness can be expected, this will hamper large-scale studies for some 
time to come. A more effective approach is to understand the important effects found 
in microscopic calculations or based on fundamental properties of nuclei and to derive a 
simple, even if phenomenological or parametrized, way to implement them in the reac- 
tion rate calculations. A good example for using microscopic (in this case, shell model 
Monte Carlo) calculations in this way is the treatment of the parity distribution in nuclei 
at low excitation energies, as discussed in Sec. 13.3.21 An example of using fundamental 
properties of nuclei is the approach to calculate nuclear level densities fSec. I3.3.T|) . 

In the NON-SMOKER code it was tried to implement reliable predictions of nuclear 
properties following the above philosophy. For instance, the nucleon optical potentials 
are calculated in a fast microscopic model [S3]- The photon transmission coefficients are 
derived using a global parametrization of GDR properties [68J, based on the hydrodynamic 
droplet model and including an energy-dependent GDR width jSHJ. Further included are 
width fluctuation corrections and isospin effects [58J. The nuclear level density will be 
discussed in the following section. Details of all nuclear properties used are given elsewhere 
PIH8] . 

3.3. Nuclear Level Density 
3.3.1. A global description 

The nuclear level density is an important ingredient in the prediction of nuclear reac- 
tion rates in astrophysics. The applicability of the statistical model of nuclear reactions 
(Hauser-Feshbach formalism) can be derived from the average level spacing and lower lim- 
its of energy and temperature for the application to calculate cross sections and reaction 
rates, respectively, can be given [pTj . 

The current version of the NON-SMOKER code uses a level density description based 
on fundamental properties of nuclei, as introduced above. The macroscopic properties can 
be derived from the Fermi-gas formalism [H]. Obviously, these will hold as long as the 
nucleons can be described by a Fermi-gas, regardless of whether they are close to or far 
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from stability. Pairing is usually accounted for by introducing a backshift, leading to the 
well-known expressions of the shifted Fermi-gas [7]: 

p(U,J,7r)=F(U,J,ir)p{U) , (1) 
with 

- 2 =^/!< u = e - s ■ 

The spin and parity dependence T is determined by the spin cut-off parameter a. Thus, 
the level density is dependent on only two parameters: the level density parameter a and 
the backshift 5, which determines the energy of the first excited state. The divergence for 
E = 5 can be avoided by either introducing an additional term depending on the nuclear 
temperature jlU] or by matching it to the constant temperature formula at low energies 

EH- 

Within this framework, the quality of level density predictions depends on the reliability 
of systematic estimates of a and 5. All current (microscopic) calculations prove that the 
backshifted Fermi-gas can account for the nuclear level density as long as a and 5 are 
chosen properly, e.g. Shell Model Monte Carlo j2H3j . combinatorial approaches and 
recurrence relations for exact level densities [7T] . 

Thus, while the overall shape of the energy dependence is given by the globally valid 
Fermi-gas, the microscopic corrections enter via 5 and mainly via a: 



a(U, Z, N) = a(A) 



l + C{Z,N)£2p- 



(3) 



where 

a(A) = aA + f3A 2/3 (4) 
and 

f(U) = l-exp(- 7 C7) . (5) 

Therefore, the parameter a contains all effects beyond a spherical droplet. Here, the 
microscopic correction C(Z,N) is thermally damped away at high excitation energies as 
suggested by |29l3Uj after inspection of hydrodynamic models. Three free parameters a, 
(3, 7 have to be fitted to nuclear data. The strength of the above approach lies in the 
fact that these microscopic properties can easily be extracted from other (microscopic or 
semi-microscopic) calculations which are available across the whole nuclear chart: nuclear 
mass models. Many different inputs can easily be implemented and compared. This is 
explained in detail in [HI]. 

Obviously, the resulting level density inherits the deficiencies of the underlying model 
concerning the microscopic properties. This is illustrated by Fig. 0] for two cases. Never- 
theless, for the majority of nuclei deviations of less than a factor of two can be obtained 
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Figure 4. Ratio of predicted to experimental level densities at the neutron separation 
energy with microscopic inputs taken from j2Z] (left) and from |49 i (right). The deficiencies 
of the underlying mass formula are propagated to the level density result. For instance, the 
well-known problems of [27| at higher mass numbers can be seen clearly. Using input from 
[4*^] improves the situation drastically. However, shell closures are still overpronounced as 
can be seen from the emerging pattern in the level density ratios. 



when using, for instance the FRDM input [IHj, which translates into a much lower uncer- 
tainty (about 30%) in the reaction rates. It has to be emphasized that the accuracy of 
this purely theoretical approach is still unsurpassed, despite recent attempts for improve- 
ment with other approaches, e.g. |22|14j (it should be noted that the average deviation 
quoted in [T3j includes results which were renormalized to experiment and therefore is 
not a measure of predictive power). 

For all details of the level density description and further implications, please refer to 
the extended paper jfiT] . 

3.3.2. Implementing parity-dependence 

So far, all theoretical, global calculations of astrophysical rates assume an equal dis- 
tribution of the state parities at all energies. It is obvious that this assumption is not 
valid at low excitation energies of a nucleus. However, a globally applicable recipe was 
lacking. We combine a formula for the energy-dependent parity distribution £Q with the 
microscopic-macroscopic nuclear level density described above |45|44j . The formula re- 
produces well the transition from low excitation energies where a single parity dominates 
to high excitations where the two densities are equal. It was tested against Monte Carlo 
shell model calculations. 

Alhassid et al. £Q have introduced a simple model for the partition function ratio 
Z-/Z+ for nuclei in the iron region using the complete pf + gg/2 shell. This model was 
combined with the shifted Fermi-gas approach to derive parity-dependent level densities. 
For details, see |45l44j . 

The parity-dependent level density was used to calculate astrophysical reaction rates 
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Figure 5. Comparison of the parity dependent reaction rates to the standard [HO], which 
assumes an equal distribution of odd and even parity states. The final nucleus is specified 
for each reaction. 



involving the three nuclides Fe, Ni, Zn in the global Hauser-Feshbach model NON- 
SMOKER. A comparison to the standard values is shown in Fig. 03 The impact on the 
rates involving the Ni and Zn nuclei is small and negligible compared to the remaining 
uncertainties in the global HF model. This is due to the fact that a sufficiently large 
number of excited states is known experimentally. Up to 20 experimental states are 
considered in the standard calculation and only above the last known state, the theoretical 
level density is in effect. However, the case is different for reactions involving 64 Fe. No 
information on experimental states is known here and therefore the full impact of the 
parity dependence can be seen. In the (n,7) case 20% difference are found. Much larger 
differences are seen in the reactions involving 64 Fe in the final particle channel. Because of 
lack of negative parities at low excitation energies, the particle emission channel becomes 
strongly enhanced in all such reactions with low or negative Q values. The (n,7) channel 
show lower sensitivity because the total transmission coefficient includes more transitions 
to states at higher excitation energy where the parity ratio is already close to unity. This 
is not true for the particle (exit) channels where preferrably states at low excitation energy 
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are populated, due to the reaction energetics. 

The case for 64 Fe shows that a large effect of the parity dependence can be expected far 
from stability where no experimental information on excited states is available and that 
it is extremely important to include such a modified level density The current approach 
is valid only for even-even nuclei in the pf + g$/2 shell. Work is in progress to extend this 
description to be able to calculate the parity distribution for a large number of nuclei far 
from stability on the proton-rich as well as neutron-rich side. 

3.4. Optical ^-nucleus potentials 

There have only been few attempts to derive global optical potentials for a-projectiles |54 
and most of them are only valid at a-energies larger than 30 MeV. Due to the high 
Coulomb barrier and nuclear structure effects defining the imaginary part of the potential 
it is difficult to obtain a global potential at astrophysical energies. Elastic a-scattering 
data can constrain the real part of the potential |42|47| and detailed analysis can also 
improve on the imaginary part |46|19j , describing the absorption into other channels than 
the elastic scattering, i.e. the Hauser-Feshbach channel. Due to the scarcity of data for 
intermediate and heavy nuclei, attempts to improve on the potential are mostly concen- 
trating on single reactions |67|19j . More global approaches suffer from the lack of data to 
confine their parameters |54|15j . 

Recently, it was tried ^Hj to find a potential for the A ~ 140 mass region by simultane- 



ously fitting data for 143 Nd(n,a) 140 Ce H>], 147 Sm(n,a) 144 Nd £H, and 144 Sm(a, 7 ) 148 Gd [ST]. 



Different parameters for the potential geometry and the energy dependence of the depth 
of the imaginary part were explored |16|55j . No significant differences were found between 
using a Brown-Rho shape g7j W{E) = W ((E - E ) 2 )/((E - E ) 2 + A 2 ) or a Fermi-type 
shape [HZ] W{E) = W Q /(1 + exp((E* - E)/a*)) of the energy dependence. For the latter, 
the parameters E* = 18.74 MeV, a* = 2.1 MeV were found, with all other parameters as 
in the previous paper (HZj- The Brown-Rho best fit was obtained with Eq = 6.35 MeV and 
A = 28.4 MeV, with the same fixed parameters V = 162 MeV, r r = 1.27 fm, a r = 0.48 
fm, Wq = 19 MeV, ry = 1.57 fm, ay = 0.6 fm. The results from the simultaneous fit of 
three reactions are shown in Figs. El and UJ 

Despite the fact that the considered targets are in the same mass region, the derived 
parameters also describe acceptably well the reaction 96 Ru(a,7) 100 Pd jSHJ. However, it 
is remarkable that even better overall agreement with all four reactions can be obtained 
when using a mass- and energy-independent potential of Saxon- Woods form for the real 
and imaginary parts (see Fig. |HJ). The real parameters are given by Vq = 162.3 MeV, 
7> = 1.27 fm, a r = 0.48 fm, the imaginary ones by Wq(E) = Wy = 25 MeV, ry = 1.4 
fm, ay = 0.52 fm. Thus, the real part is identical to the potential by Somorjai et al. jHZ] 
but without energy dependence, whereas the imaginary part is similar to the one used 
in McFadden & Satchler [HZj. Since the McFadden & Satchler parameters were derived 
from extensive elastic scattering data it seems reasonable that they are applicable to a 
wider range of targets. The Somorjai et al. parameters were derived for one reaction only 



The optical potential is parametrized as 
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Figure 6. Cross sections, reaction rates, and S-factors from a simultaneous y 2 fit of the 
Fermi-type energy- dependent a+nucleus optical potentials of three reactions (see text). 
The dashed lines are the statistical model calculation. The errors on the 147 Sm(n,«) rates 
were assumed to be 10%. 



but seem to work also for the nuclides investigated here. The new energy-independent 
potential also describes well other reactions not shown here, e.g. 70 Ge(a,7) 74 Se. 

Here, we do not show our results from fitting each reaction separately. Obviously, 
potentials fitted to single reactions can describe those - but only those - even better. 

Closer examination of Fig. |H] seems to suggest that an additional energy-dependence 
has to be introduced at very low a-energies. This can be seen mainly from the 144 Sm(a,7) 
comparison. However, since so far comparisons to other reactions were successful and the 
potentials only failed for this one case, it could also be that this is just a peculiar reaction. 
One peculiarity is the a-decay of the final nucleus 148 Gd. In fact, in the considered 
reactions, also 144 Nd is an a-emitter. As can also be seen from Figs. EHU the description 
of 147 Sm(a,n) 144 Nd is always less good than of those reactions without a- unstable final 
nucleus. This is supported by the fact that non-statistical effects have been found in 
this reaction |37|. It seems that an additional effect has to be included for a-unstable 
nuclides, either by a modified phenomenological potential or by explicitly accounting for 
the correction. In the case of 144 Sm(a,7), the effect seems to start acting at around 12 
MeV, perhaps similar to an additional barrier. 
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Figure 7. Cross sections, reaction rates, and S-factors from a simultaneous x 2 fit of the 
Brown-Rho energy-dependent a+nucleus optical potentials of three reactions (see text). 
The dashed lines are the statistical model calculation. The errors on the 147 Sm(n,a) rates 
were assumed to be 10%. 



4. Conclusions 

Stellar models have now reached a stage where the often quoted astrophysical accu- 
racy of "a factor of two" is not sufficient anymore in many cases. This poses a spe- 
cial challenge for the experimentalist as well as the theoretician. Global nuclear models 
have already been quite successful in predicting nuclear rates, especially for neutron- and 
proton-induced reactions. Despite these considerable successes close to and far off sta- 
bility, the description of certain nuclear inputs, such as optical a-potentials, still needs 
to be improved. It is also still unclear whether nuclear properties far off stability can 
be predicted with sufficiently high accuracy. Although future advances in microscopic 
theories may alleviate that problem, experimental data is clearly needed. Rare Isotope 
Accelerators will make it possible to study highly unstable nuclides but also "classical" 
nuclear physics experiments with stable or long-lived nuclei are indispensable. They are 
not only required to study specific crucial reactions in stellar evolution and nucleosyn- 
thesis but can also provide the systematics for global descriptions and shed light on the 
interaction of different reaction mechanisms. 
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